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Abstract. In this paper, by a new method we estabUsh the Weyl-type asymptotic 
formula for the counting function of biharmonic Steklov eigenvalues with Neumann 
boundary condition in a bounded domain of an n-dimensional Riemannian manifold. 



Let (M, g) be an oriented Riemannian manifold of dimension n with a positive definite 
metric tensor g, and let D C be a bounded domain with C^-smooth boundary dD. 
Assume g is a, non-negative bounded function defined on dD. Consider the following 
biharmonic Steklov eigenvalue problem with Neumann boundary condition: 



where v denotes the inward unit normal vector to dD, and is the Laplace-Beltrami 
operator given in local coordinates by 



Here \g\ := det(gij) is the determinant of the metric tensor, and g^^ are the components 
of the inverse of the metric tensor g. 

The problem was first discussed in 1968 by J. R. Kuttler and V. G. Sigillito (see 
[13] ) since it describes the deformation u of the linear elastic supported plate D under the 
action of the transversal exterior force f{x) = 0, x £ D (for example, when the weight of 
the body D is the only body force) with Neumann boundary condition I r, = (see. 



1. Introduction 



(1.1) 




in D, 
on dD, 
g^u = on dD, 




[31] or p. 32 of i33 ). 



1991 Mathematics Subject Classification. 35P20, 58J50, 58C40, 65L15 
Key words and phrases, biharmonic equation, Steklov eigenvalue, asymptotic formula, Riemannian 
manifold. 



2 



GENQIAN LIU 



It is well-known that the problem (|1.1|) has nontrivial solutions u only for a discrete 
set of A'^ ~ A|, which are called biharmonic Steklov eigenvalues with Neumann boundary 
condition. Let us enumerate the eigenvalues in increasing order: 

= A3 < A? < Ai < . . . < A3 < . . . , 

where each eigenvalue is counted as many times as its multiplicity. The corresponding 
eigenfunctions mqi, wo2, ui,U2, ■ ■ ■ ,Uk, ■ ■ ■ form a complete orthonormal basis in L^g{dD) 
(see, Proposition 3.5), where uoi{x) = 1, uq2{x) — Jj^ F{x, y)dRy on x € dD, and F{x, y) 
is Green's function in D with Neumann boundary condition. It is clear that A| can be 
characterized variationally as 

Ao — u, Ai — — ^ ^— 23 — — im —r. — ' 

}qo £' "i"^ veLip{D}nH^{D) jgjj Q-^v^ds 

Sod S^^ds = 0. Sg^ s^ugr^^ds = 

^ _ J^\AgUk\^dR _ . J^\AgV\^dR 

= p - ^— — ma,x inf ^ - r— - , k — 2, 3, 

JoD^'^kds J-c{v\vGMp{D)nHHD), "^.F Jg^g'iv^ds 

SdD 8^" dB=0, /go g^Ug2''dB = 0}, codim(^) = fc + l 

where H"'^{D) is the Sobolev space, Lip{D) is the set of Lipschitz functions on D, and 
where dR and ds are the Riemannian elements of volume and area on D and dD, respec- 
tively. 

The problem (|l.ip is also important in biharmonic analysis because the set of the 
eigenvalues for the biharmonic Steklov problem is the same as the set of eigenvalues 
of the well-known "Dirichlet to normal derivative of Laplacian" map for biharmonic 
equation (This map associates each function u defined on the boundary dD to the normal 
derivative '^^'^j"^ of A^u, where the biharmonic function u in D is uniquely determined 



In the general case the eigenvalues A^ can not be evaluated explicitly. In particular, for 



hyu\g^ and {du/diy)\g^ = 0). 

aigenvalues 

large k it is difficult to calculate them numerically. In view of the important applications, 
one is interested in finding the asymptotic formula for A| as fc — >■ oo. Let us introduce 
the counting function A{t) defined as the number of eigenvalues A| less than or equal to 
a given . Then our asymptotic problem is reformulated as the study of the asymptotic 
behavior of A{t) as t — -|-c». 

In order to better understand our problem (11.11) and its asymptotic behavior, let us 
mention the Steklov eigenvalue problem for harmonic equation 



(1.2) 



AgV = in D, 

^ + Lgv = Q on dD, 

where t is a real number. This problem was introduced by M. W. Steklov for bounded 
domains in the plane in |28J. His motivation came from physics. The function v represents 
the steady state temperature on D such that the flux on the boundary is proportional 
to the temperature. For the harmonic Steklov eigenvalue problem (|1.2p . in a special case 
in two dimensions, A. Pleijel |23j outlined an investigation of the asymptotic behavior of 
both eigenvalues ik and the eigenfunctions Vk- In 1955, L. Sandgren [55] established the 
asymptotic formula of the counting function B{t) = t^, < t}: 
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I.e., 

B{t) CJ„_1 f 

where uJn-i is the volume of the (n — l)-diniensional unit ball, ds is the Riemannian 
element of area on dD. This asymptotic behaviors is motivated by the similar one for 
the eigenvalues of the Dirichlet-Laplacian. The classical result for the Dirichlet-Laplacian 
on smooth domain D is Weyl's formula (see [3S], [33] or [S]): 

(1.4) Nd{t,D) ^ ^{vo\{D))t-/^ as r ^ 

where N{t, D) = i^ifJ.k < t} and fik is the fc-th Dirichlet eigenvalue for D. 

The study of asymptotic behavior for the biharmonic Steklov eigenvalues with Neu- 
mann boundary condition is much more difficult than that for the harmonic Steklov 
eigenvalues. It has been a tempting and challenging problem in the past 40 years. The 
main stumbling block that lies in its way is the estimates for the different kinds of Steklov 
eigenvalues corresponding to the different kinds of boundary conditions. For the simpler 
biharmonic Steklov eigenvalue problem with Dirichlet boundary condition, the author 
established the leading term asymptotic formula of the eigenvalues (see, p^). 

In this paper, for the biharmonic Steklov eigenvalues with Neumann boundary condi- 
tion, by a new method we establish the Weyl-type asymptotic formula of the counting 
function. The main results are the following: 

Theorem 1.1. Let {M,g) be an n- dimensional oriented Riemannian manifold, and 
let D <Z be a bounded domain with -smooth boundary dD. Then 

(1.5) A{t)^^^ / ^?"-lds + o(T"-l) asT^+oo, 



(^167r)"-i Jqd' 
where A(t) is defined as before. 

Corollary 1.2. Under hypothesis Theorem 1.1, if g = 1 on dD for problem il.l\) . then 

k + 2 



(1.6) Afe vIOtt — as k ^ +oo. 

\uj„-i{vol[dD)) ' 



However, when the boundary of a bounded domain is smooth, we have the following 
Weyl-type asymptotic formula with a better remainder estimate: 

Theorem 1.3. Let (A4,g) be a smooth, n-dimensional oriented Riemannian manifold, 
and let D G M be a bounded domain with smooth boundary dD. Then 

(1.7) A{r) ^ / g^^-'ds + 0{t-^) as r ^ -Hoo, 

(Vl67r)" ^ JoD 

where A{t) is defined as before. 

The proofs of our main results uses four key techniques: The first technique is the 
compact trace lemmas for the domain which is the union of a finite number of Lipschitz 
images of cubes. The second technique is to give the explicit formula for the different 
kinds of biharmonic Steklov eigenvalues and eigenfunctions in a cube of R" (by the 
method of separation variables we seek the product form of eigenfunctions, one of factors 
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is the Dirichlet eigenfuction or Neumann eigenfunction, see Section 4). Then we can use 
the well-known variational methods, which H. Weyl and R. Courant and D. Hilbert 
have employed in the case of the membrane to give the asymptotic formulas for the two 
kinds of the Steklov eigenvalues in the cube. The third technique is put the biharmonic 
Steklov problem into an abstract Hilbert space theory. That is, we first make a division 
of D into subdomains. From this division we construct two Hilbert spaces IC^ and JC^ and 
isometric mappings of into /C and /C into /C*. Those of subdomains situated at the 
boundary Tg we can map on cylinders of type treated in Section 4. In a sufficiently fine 
division of these, the variant of the g^'' and g will be small and they can be replaced by 
constants. Then, we can estimate the asymptotic behavior of the eigenvalues. By means 
of the results of Section 5, we get the asymptotic formula with leading asymptotic for 
A{t). Finally, applying Theorem 1.1 and a standard technique based on the asymptotic 
behavior of spectral function of pseudodifferential operator (see, p. 162 of [27], [11] or 
[31]). we obtain the desired result of Theorem 1.3 (with a better remainder estimate). 



2. Compact trace Lemmas 



An n-dimensional cube in R" is the set {x e M"|0 < Xk < a, k = 1, ■ ■ ■ ,n}. A set 
D C K" is said to be a Lipschitz image of a set fl C K" (see [33]) if there is a one-to-one 
map from to D defined by 

(2.1) X = ^{x'), x' en 

satisfying a Lipschitz condition 

(2.2) c-'\x' - y'\ < \^ix') - ^{y')\ < c\x' - y'\ 

for some constant c and all x' and y' in il. {\x\ — {xl + ■ ■ ■ + x^)^/^, \x'\ — {x['^ + • • • + 

x'^'y/'). 

A set D C M" is said to be star-shaped with respect to a point x^ ii x G D implies 
that the closed segment {(1 — t)x^ + tx\0 < t < 1} is contained in D. Now assume that 
Z) is a bounded domain in R" and that the closed domain D is star-shaped with regard 
to all points in an open neighborhood of a point £ D. We can assume a;° = (0, • • • ,0). 
In this section, ||x|| denotes an arbitrary norm in R" with the usual properties of a norm, 
that is 

a) \\x\\ = <^ X = 0, b) ||a; + y|| < ||a;|| + ||y|l and c) = 

where < is a real number. Then evidently (see [26]) there is a (5 > such that D is star- 
shaped with respect to all points in Bg = {a::|||a;|| < S}. Since Bs is open, it is clearly that 
X d Bs and y € D implies that all the inner points of the segment {{l — t)x-\-ty\Q < t < 1} 
belong to D. 

Lemma 2.1 (Sandgren, p 21 of [2^). // a bounded domain D C R" is star-shaped 
with respect to all points in the open cube — {x & R"| maxi<i<„ \xi\ < S}, then D is 
a Lipschitz image of the cube D' (with the side-length 2a) given by a transformation 

(2.3) X = *(a;'), xE D and x' e D' 
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satisfying the Lipschitz condition 

(2.4) c-'\\x' - y'W < \\^{x') - ^{y')\\ < c\\x' - y% 

where c = max(3a/(S, 36^/(5a) and b = niax^g/) ||a;||. 

Let / be a real- valued function defined in an open set D in R" {n > 1). For y £ D we 
call / real analytic at y if there exist e and a neighborhood U of y (all depending 
on y) such that 

for all X inU . We say / is real analytic in D, if / is real analytic at each y E D. 

From here up to Section 5, let M be an n-dimensional Riemannian manifold with real 
analytic metric tensor g. We say that D is a Lipschitz image of a cube if it is contained 
in some coordinate neighborhood U and its image Di in given by the coordinates of 
[/ is a Lipschitz image (see, previous definition) of a closed cube in R". 

A subset F of {A4,g) is said to be an (n — l)-dimensional smooth surface if F is 
nonempty and if for every point x in F, there is a smooth diffeomorphism of the open unit 
ball B{0, 1) in R" onto an open neighborhood U of a; such that B{0, l)n{x € M"|a;„ = 0} 
maps onto U CiT. 

Let D together with its boundary be transformed pointwise into the domain D' to- 
gether with its boundary by equations of the form 

(2.5) x'i ^ Xi + fi{xi,- ■ ■ ,Xn), 2 = l,2,---,n. 

where the functions fi and their first order derivatives are Lipschitz continuous through- 
out the domain, and they are less in absolute value than a small positive number e. Then 
we say that the domain D is approximated by the domain D' with the degree of accuracy 
e. 

It is well-known (see, for example, p. 133 of [10 or p. 24 of [26]) that every element u 
in Lip(D) has partial derivatives du/dxk, k ^ I, ■ ■ ■ , n, which are defined a.e. in D and 
belong to L"°{D). In particular, Lip{D) C H^{D). 

A subset ^ of L^{dD) is said to be precompact if any infinite sequence {uk} of elements 
of contains a Cauchy subsequence {uk'}, i.e., one for which 



/ (life' — M;')^(is — )• 0, when k'J'^oo. 

JdD 



Lemma 2.2. Let D C {M,g) be a Lipschitz image of a cube and let g be a non- 
negative function in L°^{dD) such that g^ds > 0. Assume 9Jl is a set of functions u 
m N{D) = {u\ue Lip{D) nH^{D), = on dD} for which 

(2.6) / \Au\^dR+{ [ Q^udsf 

J D JdD 

is uniformly bounded. Then the set : m G 371} is precompact in L^{dD). 

Proof. It follows from Theorem 2.3 of [26 that there exists a constant Ci > 0, only 
depending on D and g, such that for all u € Lip{D), 



(2.7) 



/ u^dR<Ci I \S/gufdR+{f g^udsf 

Jd Ud JdD 
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where 



ID J4>{D) 

and (j){D) is the coordinate image of D. Put 

v2 



J^\Agu\^dR+{jQ^ g^udsY 



(2.8) Ai{D) = inf 

ueN(D) f^\Vgu\^dR+ ( L^g^udsY 

In order to prove the existence of a mininiizer to (|2.8p . consider a minimizing sequence 
Um in the set iV(Z?), i.e., 

/ \AgUrn\'^dR + / ds)^ Ai(_D) as m — !• +00 

Jd JdD 

with |VgU„ip(ix + ( Jgjy g^u„i ds)^ — 1- Thus, there is a constant C2 > such that 

(2.9) \\AgU^\\L-2(D)<C2, / iVgU^pdx < C2, (/ ^^V.ds)'<C2 

J D JdD 

for all m > 1. It follows from the a priori estimate for elliptic equations (see, for example. 
Proposition 7.2 of p. 345 in [501) that there exists a constant C3 > depending only on 
n, D such that 

(2.10) \\Ur,i\\m[D) < CiiWAgUmllL-^iD) + 1 1 || ffi (D) ) • 

From this, (PTT]) and we have that 

||Um||H2(£,) < C4 for all TO. 

By the Banach-Alaoglu theorem we can then extract a subsequence, which we still call 
{um}, converging weakly in H^{D) to a limit m, and strongly converging to u in L^{D). 
Since the functional \Agu\'^dR is lower semicontinuous in the weak H^{D) topology, 
we have 



/ \AgufdR < lim / \AgU„,\^dR, 

J D m— >oo J D 



Since u,„ — >■ u weakly in H (D), we get that u„i — > u strongly in H'^(Di) for any 
< r < 2. Note that %^L^ = 0. It follows that giil ^ = 0. Therefore u e M is a. 

ov IdD Of IdD 

minimizer. 

We claim that Ai{D) > 0. Suppose by contradiction that Ai{D) — 0. Then 



(2.11) 

and 



AgU = in Z?, 
^ = on dD 

av 



(2.12) / g'^uds^Q. 

JdD 



The boundary value problem (12. lip implies that u = constant in D. By J^^ g^ ds > Q 
and ()2.12p . we then get u = in 13. This contradicts the fact that jVgUpdi? + 

( /a_D ^''^^ rfs)^ — 1, and the claim is proved. 
By (|2.8p . we obtain that for every u e N{D), 

(2.13) f\S/gufdR+{f g\dsy <-^( f \AgufdR+{ [ g^dsf] . 
Jd JdD ^iW) \Jd JdD J 
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It follows from ([TT]) and ([2?T3| that is a bounded set in H^{D). Since D C {M,g) is 
a Lipschitz image of a cube, it follows from |^ (see also, Chs V, VI of [6J that the set 
{u\qj^ : u £ EOT} is precompact in L^{dD). □ 

Lemma 2.3. Let {M,g) be a real analytic Riemannian manifold, and let D C {M,g) 
be a Lipschitz image of a cube. Assume ttiat Ti is a portion of dD and ttiat Tq is an 
{n — 1)- dimensional C'^'^ -smooth surface in dD satisfying Tq <Z<Z dD — Fi. Suppose that 
g is a non-negative function defined on dD and assume that € ( respectively, &) is a set 
of functions u in K'^{D) = {u\u G Lip{D) D {D) , = on Ti, u = ^ = on Tq} 
(respectively, K{D) = {u\u € Lip{D)r)H'^{D), ^ = on TiUTq, and u = on dD- 
Ti}) for which 

(2.14) / \Agu\^dR+( [ g^uds 



D \JdD 



is uniformly bounded. Then the set : u G £} (respectively, : u G &}) is 

precompact in L^{dD). 

Proof. We only prove the case in K'^{D) because the method is similar for the case in 
K{D). It follows from Theorem 2.3 of 26 that there exists a constant C > 0, depending 
only on D and g, such that for every u G Lip{D), 

2" 



|V„uPdi?+ / g^uds 



dD 



(2.15) / u^dR < C 

Jd 

Let 

. , fn\^nv\^dR+ (L^g^vds? 
2.16 Ap,(D) = inf Jd\ 9 I VJaD£ J 

r ,v ... , J^\^,v\^dR+{J,^Pvdsy 

In order to prove the existence of a minimizer to (j2.16l) , consider a minimizing sequence 
in K'^iD), i.e., 

J \AgVjn\^dR+ g^Vmds^ Ag{D) as rrn- +oo 

with \\7 gVml'^dR + {jgjj g^Vm ds)^ = 1. Thus, there is a constant C > such that 
(2.17) ||AgW,„|||2(B) + ^^e^'^ds^ <C, 

l|Vs«m||i2(£,) + g^Vmds^ <C for ah m. 

Let {Di} be a sequence of Lipschitz domains such that Di C D2 C ■ ■ ■ C Di C ■ ■ ■ CC 
{DUTiU To), W^iDi = D, and Fi U Fq C dDi for all /. It follows from the a priori 
estimate for elliptic equations (see, for example, the proof of Proposition 7.2 of [30]) that 
there exists a constant C; > depending only on n, Di,D, Fi and Fq such that 

(2-18) hrnWmiDi) < C'liW^gVmWmD) + 1 1 1 1 ffi (D) ) • 

From this, (P?T5|) and (PTT)) . we have that 
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SO that 



< 



where C'/ is a constant. For each /, by the Banach-Alaoghi theorem we can extract a 
subsequence {vi^m}m=i of {vm}, which converges weakly in H^{Di) to a limit u, and 
strongly converges to u in L^{Di). We may assume that {vi+i^m} is a subsequence of 
{vi.m} for every I. Then, the diagonal sequence {vi,i} converges weakly in to u, 
and strongly converges to u in L^, in every compact subset E oi D. It is obvious that 
ll^ff^lli^j^)) + {Jgjj Q^udsj — 1. Since the functional \/Sgu\^dR is lower semicontin- 
uous in the weak H^{Di) topology, we have 

/ |AgUpdi?< lim [ \AgVkM^dR, 

J Di fe— >cxD J Di 

I \Agu\^dR = lim / lAgul^di? < lim ( lim / \AgVk.k\^dR 

Jd 1^°°Jdi l^°°\k^ooJDi 

lim ( lim / lAgVk.kl^dR) ^ \jm \AgVk,k\^dR. 

In addition, for each fixed /, since Vk,k u weakly in H^{Di), we get that Vk.k — >■ u 
strongly in H''{Di) for any < r < 2. Note that %7^|p^ = and Wfe,fe|p^ = ^^l-p^ = 0. 
It follows that = and uL = =0. Therefore u e K'^(D) is a minimizer. 

Cl^ \ L I lip ' -L 

We claim that Ag{D) > 0. Suppose by contradiction that 

lD\^9^\^dR+{Jg^pudsy 

It follows that AgU = in Z?. Since the coefhcients of the Laplacian are real analytic 
in D, and since Fq is a C'^'^-smooth surface, we find with the aid of the regularity for 
elliptic equations (see. Theorem A of [12], [H] or |T]) that u is C^'"^ -smooth up to the 
partial boundary Fq. Note that u = ^ = on Fq. Applying Holmgren's uniqueness 
theorem (see. Corollary 5 of p. 39 in [53]) for the real analytic elliptic equation AgU = 

in D, we get that m = in D. This contradicts the fact |VgUpc?i?+ {Jqj^ Q^u ds) = 1, 
and the claim is proved. Therefore we have that 

(2.19) [ \Vgu\^dR+([ g^uds^ 

JD \JaD 



1 

< — 
- A, 



\Agu\^dR+ ( / g^uds 



D \JdD 



for u e K'^{D). 



According to the assumption, there is a constant C" such that 

II'^9'"IIl2(u) + g^uds^ <C" forallwee, 

and hence 

(2.20) \\^g'u\\h{D)+ (^J fuds^ <C" for all wee. 

Combining this and ()2.15p . we have 

Hull HI (D) < C, for aU u G €, 
which implies that : m G £} is precompact in L^[dD). □ 
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Corollary 2.4. Lemma 2.3 is still true if we exchange {2.14^ for 
(2.21) / lAgUpdi?. 

JD 

Proof. Let q G L°°{dD) be the characteristic functions of Tq (respectively, F — Fi) in 
the case of K'^{D) (respectively, K{D)). Then Jp g^uds = for all u in dJl, and hence 
according to (I2.2ip . the condition (|2.14l) is satisfies such that Lemma 2.3 can be applied. 
□ 

Let Di, D2 C (M, g) be Lipschitz images of cubes with boundaries dDi and dD2- We 
say an open domain D with boundary dD and closure D = D U dD is composed of Di 
and D2 if 

a) dDi n dD2 has positive measure 

b) D = D1UD2 

c) dD c dDi U dD2 

d) dD has positive measure. 

The requirement d) excludes for instance the possibility that M is a sphere and D = M 
the union of two hemispheres Di and D2 (see, p. 27 of }26l). 

By a finite number of domains, each of which is a Lipschitz image of a cube, we can 
obtain more domains according to the above method. Denoted by J- all such domains. 
Completely similar to the proofs of Lemmas 2.2, 2.3, we find that the compact trace 
Lemmas 2.2, 2.3 are also true for each domain in class J^. 



3. Some completely continuous transformations and their eigenvalues 

Let {Ai , g) be an n-dimensional real analytic Riemannian manifold and let D C M be 
a bounded domain with boundary F. Suppose that D is of the type defined in Section 2 
(i.e., D G T) so that the compact trace lemmas 2.2, 2.3 are true. Let ghe a, non-negative 
bounded function defined on F or only on a portion Tg of F (measure F^ = /p ds > 0) 
and assume that Jp g^ds > 0. In the case Fg ^ F we let Fq be a C^''^-smooth (n — 1)- 

dimensional surface in F — Fg. 

If Fg 7^ F (measure F — F^ > 0), we denote 

K{D) = {u\u e Lip{D) n H'^{D), = on F^, U Fq, and u = on F - F J, 
K'^{D) = {u\u e Lip{D) n H^{D), = on Fg, and u = |^ = on Fq}. 

If Fg = F, we denote 

- du f f 

N{D) = {u\u e Lip{D) n H'^{D), — = on F, / g^uuoi ds ^ and / g^uu^^ ds = 0}, 

where uqi — 1 and uq2 are two eigenfunctions corresponding to the Steklov eigenvalue 
A = (see, Section 1). We shall also use the notation 

{u,v)* ^ [ {Agu){Agv)dR, u,v&K{D) or L'^iD) or N{D). 
Jd 
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The bilinear functional {u,v)* can be used as an inner product in each of the spaces 
K{D), K'^{D) and N{D). In fact, {u,v)* is a positive, symmetric, bilinear functional. 
In addition, if {u,u)* = 0, then Agtt = in Z?. For u e K{D) or N{D), by applying 
Green's formula, we have 

= / u{Agu)dR^- [ \Vgu\^dR~ [ u^ds = - [ \Vgu\'^dR, 

JD JD JdD JD 

which implies u = constant 'm D. In the case u £ K{D), in view oi u — onF — Fg, 
we get that m = in D; In the case u e N{D), since g^uds = with Jp ds > 0, 
we obtain u = in Z?; In the case u £ K'^{D), since Fq is a C^''^-smooth surface and 
u = ^ = on Fq, we find by Au = in D and Holmgren's uniqueness theorem (see. 
Corollary 5 of p. 39 in M) that u = in D. Closing K{D), K'^{D) and N{D) with 
respect to the norm ||u||* = {u, u)*, we get the Hilbert spaces {)C, \\ ■ \\*), {JC^, \\ ■ \\*) 
and (A/", || • ||*), respectively. 

Next, we consider two linear functionals 

[u, v] — g^uv ds 

and 

(3.1) {u,v) = {u,vy + [u,v], 

where u,v E K(D) or u,v E K'^{D) or u,v E N{D). It is clear that {u,v) is an inner 
product in each of the spaces K{D), K'^(D) and N{D). 

Lemma 3.1. The norm 
and 

\\u\\ = y/(u~u} 
are equivalent in K{D), K'^{D) and N{D). 

Proof. We have to show that there exist positive constants Ci and C2 such that 

Ci||u||* < ||u|| < Csllwir for aU u in K{D) or K'^{D) or N{D). 

Obviously, ||w||* < ||m|| for all u in each of three spaces. Let us first consider the case 
u E N{D). It suffices to show that ||u|| is bounded when u belongs to the set 

m = {u\u E N{D), \\u\\* < 1}. 

It follows from Lemma 2.2 that 9Hr := {w|p : u E OJl} is precompact in I/^(F). This 
implies that there exists a constant C > such that Jp u'^ds < C for all u E 9Jl. Therefore, 
[u,u] = j^g^v?ds is bounded in OT, and so is — {u,u)* + [u,u]. 
Next, we consider the case for u E K{D) and still denote 

m={u\uE K{D), \\u\\* < 1}. 

By taking 

. _ r 1 for X e F - Fg, 
^ ~ \ for X e Tg, 



ASYMPTOTIC FORMULA FOR STEKLOV EIGENVALUES 



11 



we have 

(3.2) / \Au\^dR+{ I Q^udsf <C for every u e »l. 

JD Jr 

It follows from Lemma 2.3 that 3Jtr is precompact in L^(r). In particular, /p ds is 
bounded on 9Jl and hence also [u,u] = /p g^u^ ds and = + [u,u]. 

Similarly, by taking g to be the characteristic function of Fq and by applying Lemma 
2.3, we can prove the corresponding results for the space K'^{D). □ 

From Lemma 3.1, it follows that 



\[u,u] 



g'^u^ds 



<C{u,u)* for all u in K{D) or K'^{D) or N{D). 



Therefore, [u, v] is a bounded, symmetric, bilinear functional in {K(D), {■,■)*), {K'^{D), {■,■)*) 
and iN{D), {■,■)*). Since it is densely defined in (/C,(-, •)*), (JC^, (•, •)*) and (A/", (•, •)*), 
respectively, it can immediately be extended to (/C, (•,•)*), {IC^, {■,■)*) and (A/", (•, •)*). 
We still use [u, v] to express the extended functional. Then there is a bounded linear 

transformation of (/C, (•, •)*) into (/C, (•,•)*) (respectively, G^i of (/C'', (•, •)*) into 

(/C^ (•, •)*), G^J}^ of (A/", (•, ■}*) into (A/", (•, ■}*)) such that 

(3.3) [u,v] = {Gl^\,vy for all u and w in /C 
(respectively, 

(3.4) = (G^iu,t;)* for aU u and v in /C*, 

(3.5) [w,!;] = (Gj^^u,w)* for all u and v in TV). 

Lemma 3.2. The transformations G\^\ G^] anrf G^j^'^ are self-adjoint and compact. 

Proof. Since [u,v] is symmetric, we immediately get that the transformation g\^\ 

G^i and Gy^'* are all self-adjoint. For the compactness, we only discuss the case for the 

transformation G^'. It suffices to show (see, p. 204 of 25J): From every sequence {um} 
in K(D) which is bounded 

(3.6) lI'Mmll'^ < constant, m = 1, 2, 3, • • • , 
we can pick out a subsequence {um'} such that 

(3.7) (gW( 

Applying Lemmas 2.3, 3.1 with the aid of (13.61) . we find that the sequence {umL } is 
precompact in L^(Fg), so that there is a subsequence such that 



/ {um' — ui')'^ ds —i' as m',Z'— >oo. 



Therefore 



[u,n' ~ Ul' , Um' -Ul']= ("m' - Ul'f ds BS m.\l' 

which implies (|3.7p . This proves the compactness of G^''. □ 
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Except for the transformations G^'' , G^] and G')^ , we need introduce corresponding 
transformations Gk;, G]Qd and G^v by the inner product (•, •). Since 

(3.8) < [u,v] < {u,v) for all u in K{D) or K'^{D) or N{D), 

there is a bounded linear self-adjoint transformation G/c of (/C, (•, •)) (respectively, Gj^d 
of (/C^, (•, •)), Gj^ of (A/", (•, •))) such that 

(3.9) [u,v] = {Gjcu,v) for all u and u in /C 
(respectively, 

(3.10) = {Gx,dU,v) for all u and v in /C^, 

(3.11) = {G_\fu,v) for all u and u in A/"). 

Lemma 3.3. The transformations Gjc, Gf^d and Gj\f are positive and compact. 

Proof. From > for any it e /C or JC^ or A/", we immediately know that Gic^ 

GjQd and G^ are positive. The proof of the compactness is completely similar to that of 
Lemma 3.2. □ 

It follows from Lemma 3.3 that Gfc (respectively, G^d, Gjv") has only non-negative 
eigenvalues and that the positive eigenvalues form an enumerable sequence {y^K.} (respec- 
tively, {/i^;''}? {^-m}) with as the only limit point. 

Theorem 3.4. The transformations G^' and G/c (respectively, G^] and G/cd, G^j^ 
and Gj\f ) have the same eigenfunctions. If and ^jc ( respectively, /ij^^ and ^^^d , fi^ 
and ^j\f ) are eigenvalues corresponding to the same eigenfunction we have 

(3.12) ^i,c = ^^'^ 



( respectively, 



(3.13) ^i^d - 



1 + m;^. ' 



Proof. We only prove the case for the Gjc (the arguments are similar for Gj^d and G_\f). 
Since G^'' is positive, we can easily conclude that the inverse (1 -I- G^'')"^ exists and is 
a bounded self-adjoint transformation. By virtue of p.3p . p.9p and (|3.1|) . we have 

(3.15) {G^^^u,vy = [u,v] = {G,cu,v) 

= {G,cu,vr + {G'j^^G,cu,vr 

= {GKU,vy + {GkG\^K,v)\ iu,velC). 

It follows that 

(3.16) G^ = gW(1 + G«)-\ 
from which the desired result follows immediately. □ 
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Proposition 3.5. Letu andv be two eigenfunctions in (/C, (•, •)) (respectively, {IC^, (•, •)), 
{Af, {■,■))) of the transformation Gk: (respectively, Gi^d, Gj\f) at least one of which cor- 
responds to a non-vanishing eigenvalue. Then u and v are orthogonal if and only if the 
u\„ and U p are orthogonal in L^(r„), that is, 

(3.17) [u,v] = / g^uvds = 0. 



Proof. Without loss of generality, we suppose that u is the eigenfunction corresponding 
to the eigenvalue 7^ 0. Then 

[u,v] = {Gk.u,v) = ti{u,v), 

which implies the desired result. □ 

We can now prove 

Theorem 3.6. Let D C {A4,g) be a bounded domain with piecewise smooth boundary 
r, and let D £ J- . If u is an eigenfunction of the transformations or Cj,^' with 

eigenvalue /i* 7^ 0, then u has derivatives of any order in D and is such that 

A^u ^0 in D, 
. TT ^0 on r^UTa, u^OonT-T^, 
^•^■^^l ^ AgU = Q on T - {Tg\JTo), 

- X^ff'u = on r^, with A3 = ^. 



Proof. Since u is an eigenfunction of (i.e., G^^u = fi^u), we have that ^ ' 



and M L p =0, and that 



/ Q^uvds = ^i* i {Agu){Agv)dR for all w e /C(L'). 

JTg JD 

Applying Green's formula (see, p. 114-120 of [15], [3]) to the right-hand side of the above 
equation, we obtain that 

^ / g'uvds^ I {Alu)vdR- f {Agu)^ds+ f ^^^^vds 

for all V e 1C{D), where G H-i{T) (see [15]). From ^^^^ = and v\^_^ = 0, 

we get 

(3.19) / {Alu)vdR- f {Agu)^ds+ f 

Jd Jr-(r,uro) Jr, 

for ah V e K.{D). By taking aU v e C^{D), we have A^u = in D. It follows from 
the interior regularity of elliptic equations that u € C°°(D). Noticing that uL and 
fi^lr-(r uFo) throughout space L'^{Tg) and L^(r — (Tg U Fq)), respectively, when v 



djAgu) 1 3 1 , n 

OV pL* 



runs throughout space K{D), we see that 

di' p* 

Therefore, p.lSp holds. In a similar way, we can prove the desired result for Gj^f. □ 



Agu = on r - (Tg U To), and \^ ' -g^u on Tg. 
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Theorem 3.7. Let {A4,g) he a real analytic Riemannian manifold, and let D C {A4,g) 
be a bounded domain with piecewise smooth boundary T. Let D ^ !F . Assume that Fq 
is a C^''^ -smooth (n — 1)- dimensional surface inT — Tg. If u is an eigenfunction of the 

transformations G^] with eigenvalue /x* ^ 0, then u has derivatives of any order in D 
and is such that 

in D, 



(3.20) 



Alu = 

on r„ 



u=^ = on To, 
Agu = 0, %^ = on r - (Fg U To), 



dv 



- £1^ M = on Tg, with = 4r- 



Proof. IfG^^M 
that 

(3.21) 



IJL*u, then we have that ^ = on Fg and 



du 



on Fq, and 



/ Q^uvds = ^i* i {Agu){Agv)dR for all v e /C'^(D), 



By using Green's formula and noticing that |^ I „ 

OV I i „ 



and = 1 = 0, we get that 



I To 



(3.22) 



iAlu)vdR- 



-I 



dv 
dv 



— Q^u] vds+ [ 
f^* ) Jr 



r-(r^uro) 



dv 



vds 



r-(reuro) 



'Agu)^ ds = for all v G K'^{D), 



where e iJ^i (F - (F^ UFq)). By taking all u € C^{D), we obtain that A^u = 

in D. Note that ^|p_^p ^-^^^ and i'|p_pjj run throughout the spaces L'^{T — {Tg U Fq)) 
and L^(r — Fo), respectively, when v runs throughout the space K'^{D). Thus we have 

d{AgU) 



AgU 



and 



dv 



on F- (FgUFo), 



and 



dv 



Q^u = on Tp. 



a 



Theorem 3.8. Let {M,g), D, Tg and To be as in Theorem 3.7. Assume that 7^ and 
k;| are the k-th Steklov eigenvalues of the following problems: 



(3.23) 



and 



A^u = in D, 



U - 

du 



dv 



1^ = on To, 
and =0 on F - (Fg U Fq), 



- <r £> M = on Fg 



Alu = in D, 



du 
dv 



on Tg 

du 



M=§^ = on Fo, 
AgU = and =0 on F - (Fg U Fq), 



on F 
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respectively. Then <;j^ < k| for all k > 1. 

Proof. For < a < 1, let Wfe = Uk{a, x) be the normalized eigenfunction corresponding 
to the k-th Steklov eigenvalue Afe for the following problem: 

' A^Uk = in D, 
^ = on r„, 
"fe = It = on To, 



(3.24) 



aAgUk + (1 - = and ^^^^ 

2i^-Xg^u, = on T,. 



on r- (r^UTo), 



It is easy to verify (cf, p. 410 or Theorem 9 of p. 419 in [5j) that the k-ih. Steklov 
eigenvalue \k = Afc(a) is continuous on the closed interval [0, 1] and differentiable in the 
open interval (0, 1), and that Mfc(a, x) is also differentiable with respect to a in (0, 1) (cf. 
[8]). We will denote by ' the derivative with respect to a. Then 



(3.25^ 



=0 in D, 



on r. 



dv 

"fe = "air = on To, 



duk 
dv 



(l-")^=Oand 



AgUk + aAgU'f. 

I i%<I-A',^tf"-A,3M=0onr 



9(A 



on r - (Fg U To) 



Multiplying p.25p by Ufc, integrating the product over D, and then applying Green's 
formula, we get that for < a < 1 







Jd Jd JdD Cl^ 



, d{AgUk) , 

Ui. ds — 

OD dv 



,d{Auk) 
Ui. — ds — 



^i^ds+l {Agu',r-^ds 



dD 



dv 



dv 



r-Cr^uro) 



{AgUk)^ds 
ov 



djAu',) ^ 
Uk — ds 



Xg^Uk u'f. ds 



{AgUk)^ds 
r-(reuro) 

1~ a duk \ du'f. 



r-(reuro) 
+ / (^—X'g^Uk — Xg'^u'i. ) Uk ds 



dv I dv 



ds 



r-CPgUFo) 



1 

a 



1 duk 



1 — a du'i. \ duk 



-- AgUk + - ^ - ^-^ l^\l^ds 



-A' 

-A' 



/ g'^uf. ds 4 


f \ 




Jr-(rgUFo) L 


/ g^ulds + 


/•-(FgUPo) ( 


1^, 





a dv a dv 

1 — a\ duk 
) dv 

-^-^Xds, 



dv 

}_ duk 
a dv 



duk 
dv 



ds 
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I.e., 



A',(a) = ^r^(r.uro) V. 9. ; ^ ^ for all < a < 1. 

This implies that Xk is increasing with respect to a in (0,1). Note that if we change 
the a from to 1, each individual Steklov eigenvalue Afc increase monotonically form the 
value <,k which is the fc-th Steklov eigenvalue of (|3.20l) to the value Kk which is the k-th 
Steklov eigenvalue p. 241) . Thus, we have that i^k < i^k for all k. □ 

Conversely, we can show that a sufhciently smooth function satisfying (13.18^ (respec- 
tively, p.20p ) is an eigenfunction of or cj^'' (respectively, G^^i). 

Proposition 3.9. Let D be bounded domain with piecewise smooth boundary, and let 
D E T ■ Assume that u belongs to C^{D) and let A > 0. 

a) If Tg T and u satisfies h3.18]) . then u ^ K. and u is an eigenfunction of G^'' 
with the eigenvalue /i* — X^^ , 

(3.26) G^*^u = X-^u 

is 

with the eigenvalue fj,* — X^^, 

(3.27) G^*lu = X-^u. 

c) If^g = r and u satisfies \3.18\) . then u E M and u is an eigenfunction of G^j^^ 
with the eigenvalue n* = X~^, 

(3.28) Gj^'u = X-^u. 



b) IfTgy^T and u satisfies i3.2(J\) . then u G IC^ and u is an eigenfunction of G^^^l 



Proof. i) Tg 7^ r. We claim that there is no eigenvalue A'^ = 0. Suppose by 
contradiction that there is a function u in C^{D) satisfying 

r A^u in D, 
(3.29) i f^^OonTeUro, u = on T-Tg 

[ Agu = on r - (r^ U To), and = on Tg. 

Multiplying the above equation by u, integrating the result over D, and then using 
Green's formula, we derive 

= / u{Alu)dR^ I \Agu\^dR- I u^^^^ds 
Jd Jd Jr ov 

+ / {AgU)^dS = / \AgU\^dR. 

Jr c*^ Jd 
This implies that AgU = in Z), so that 

= / u{Agu)dR^~ [ \Vu\^dR- [ u^ds^- [ \Vu\^dR. 
Jd Jd JdD c/^ Jd 

That is, u = constant in D. Since u = on F — Tg, we get that u — in D. The claim 
is proved. 
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In view of assumptions, we see that m G /C. By p.lSp and Green's formula, it follows 
that for an arbitrary v G K{D) 



= A 

= 



-3 f diAgu) 



dv 



vds = 



dv 



V ds 



dv 



dv 



{Agu) — ds+ / {Agu){Agv)dR - / w(A>)di? 



D 



D 



{Agu){Agv)dR = X^-\u,v)* 



Therefore, 



{gI^^u - X-^u, vy ^ for all v £ K{D), 



which implies (j3.26p . By a similar way, we can prove b). 

ii) Tg = r. In this case, for the eigenvalue A'^ = 0, the problem (I3.18P has the 
solutions 1*01 — constant and 1/02(2;) = SjjF{x,y)dRy in Z?, here F(x,y) is Green's 
function with Neumann boundary condition (see. Section 1). These solutions do not 
belong to N{D). If, however, m is a solution with eigenvalue A^^ > then u £ Af. Indeed, 
by Green's formula we get 



dD dv 



ds = / A^u dR = Q 

D 



and hence form p.lSp we obtain 



dD 



g'^uuQids = I Q'^uds — 0. 



dD 



In addition, from (13.18^ we get 



Ag{Agu) = in D, 

^ x^gSu on dD. 



so that 



Agu{x) = / F{x,y) 



d{AgU) 
dVy 



IdD ^"y 
Combining this and Green's formula, we have 



dsy = X^ F{x,y)g^{y)u{y)dsy 



dD 







du 
—ds 
dD ov 



AgudR = A'^ 



D \JdD 



F{x,y)g-\y)u{y)dsy dR.. 



A' / g'iyHy) [ I F{x,y)dR, ] dsy 

JdD 



dD 



D 



X / g {x)u{x) [ / F{x,y)dRy \ dsx — X" / g {x)u{x) uo2{x) ds. 



dD 



i.e., fg^ g^UQ2uds — 0, so that u £ J\f. Proceeding as in a), we can prove that (|3.28p 
holds. □ 
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Remark 3.10. Each of transformations Gj^, G]^d and corresponds to a biharmonic 
Steklov problem given by the quadratic forms 

{u,uy ^ / lAgWpdi? 

J D 

and 

and the function classes of /C*, tC^* and A/"*, respectively. The eigenvalues of these 
biharmonic Stcklov problems are given by 

(3.30) A3 = 1/^^, fc = l,2,3,---. 

Since is the only limit point of the only possible limit points of A| are +cx). 



4. Biharmonic Steklov eigenvalues on a rectangular parallelepiped 



Let D ~ {x ^ R"|0 < < Zi, i = 1, • • • , n} with boundary F, and let Tg — {x ^ 
K"|0 < Xi < li when i < n, Xn — 0}. Let F'" = {x G R"|0 < Xi < k when i < n, a;„ = 
Our first purpose, in this section, is to discuss the biharmonic Steklov eigenvalue 
problem on rt-dimensional rectangular parallelepiped D: 



( A^u = 

1^-0 



(4.1) 



u 



du 



u = Au^O 

djAu) 
du 



in D, 
on Fg, 
on F'", 

on F- (FgUF'"), 



\ g u = on Fg, g — constant > on Fg. 

We consider nonzero product solution of (|4.ip of the form: 

u = X{xi,--- ,x„_i) r(a;„), 

where X{xi, ■ ■ ■ ,a;„_i) is a function of variables xi, ■ ■ ■ ,Xn~i and Y(xn) is a function 
of Xn alone. Since 

Au = {An-lX{xi,--- ,Xn-l))YiXn)+2VX{xi,--- ,Xn-l)-VYiXn) 
+ {X{X1,--- ,Xn-l))Y"{Xn) = {An^iX{xi,--- ,Xn^l))Y(Xn) 

+ {Xixi,--- ,x„_i))y"(a;„) 
A^u = (A2_iX(a;i, • • • , a;„_i))y(x„) + 2(A„_iX(a;i, • • • , a;„_i))y"(a;„) 



and 



where 



+ {X{xi,--- ,Xn^i))Y""{Xn), 
A„_iX(xi, 



, _ ;p d'x 

i—1 
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we find by A^u that 

{A^X{xi,--- ,Xn-i))Y{xn) + 2iAXixi,--- , x„_i))r"(a;„) 
+ {X{xu--- ,x„_i))y""(x„) = 0, 

so that 

X{xi,- ■ ■ ,Xn-l) X{xi,- ■ ■ ,Xn~l) Y{Xn) Y{Xn) 

Differentiating (14.21) with respect to Xn, we obtain that 



AX{xi, ■ ■ ■ ,Xn-l) 



X{xi 



Y'\x„) 



YiXn) 



Y""{Xn) 



YiXn) 



0. 



The above equation holds if and only if 
AX{xi, ■ ■ ■ ,x„_i) 



(4.3) 



X{xi, ■ ■ ■ ,Xn-l) 



\Y""(x„)] 


/ 


[ Y{x„) \ 




2 




/ 









where r]"^ is a constant. Therefore, we have that 

(4.4) AX{X1, ■■■ , Xn^l) + V^X{X1, ■■■ , Xn-l) = 

and 



r""(x„) 



Y"{x,,) 



Y{Xn) 



= 0. 



YiXn) 

From (14. 4p . we get 

(4.5) A^X = -rj^AX ^ ifX. 
Substituting this into (|4.2p . we obtain the following equation 

(4.6) Y"'\x^) - 2rfY"{xn) + V^Y{xn) = 0. 

It is easy to verify that the general solutions of (|4.6p have the form: 

(4.7) Y{xn) = A cosh 7/x„ + B sinh riXn + Cxn cosh rjXn + Dxn sinh rjXr, 
By setting y(0) = 1, r(/„) = 0, r'(0) = 0, y'(/„) = 0, we get 

(sinh r//„) (cosh ry/„) + ry/„ 



(4.8) 



Y{x^) 



cosh "qXr, 



sinh rjln — ffln 



sinh r]Xr, 



7]{smhr]l„){coshr]ln) +rflr. 



sinh rjln — ij'^l^ 



cosh T]Xn 



rj sinh -ql^ 



sinh rjln — V^^n 



sinh rjXr, 



It is well-known that for the Dirichlet eigenvalue problem 

, , ( AX{xi,--- ,Xn-i)+ri'^Xixi,--- ,Xn-i) ^0 in D, 

\ u = Q on d{{xi, ■ ■ ■ ,x„_i)|0 < Xi < h, i ^ 1, ■ ■ ■ ,n - 1}, 

there exist the eigenfunctions 



(4.10) A (xi, • • • ,x„_i) = c ( sm — xi 



Sm — ; Xn-l 



In- 



20 



GENQIAN LIU 



which correspond to the eigenvalues 

n-l 



n-l / V 2 

rriiTT 



i=l 



where to; = 1, 2, 3, • • • 



Therefore, 
(4.11) 



miTT 
c I sm — : — Xi 



sm — Xn-l 



f (sinhry/n)(cosh?//„) + 
V sinh^ rjln - v'^l^ 
r]{smhr]ln){cos'hr]ln) + rfln 



sinh rjln — V^^r, 



[cosh rjXn 
sinh rjXn 

coshrjXn 



/ 1] sinh^ rjl„ 
Vsinh^ rjln — 'rfln 



sinh rjXr, 



Since 

we obtain 
djAu) 
dv 

so that 
with 



r"'(0) = 2r/3 ^ (sinh ?7;„) (cosh r//„) +7?;^ 



sinh rjln - rfl'^^ 



,„^o = (A„_iX(a;i,--- ,x„-i))r'(0) + (X(xi 



2Tf 



/ (sinh jy^n) (cosh ?7^„) + r]lr, 



sinh ijln — rj^l^ 
diAu) 



X{xi, 



,Xn^l))Y"'iO) 
, Xn—l), 



dv 



X'^g'^u = on r. 



^3 _ i-q^l^ / (sinh ?7?n) (cosh ?jl„) + ^„?? 

Our second purpose is to discuss the biharmonic Steklov eigenvalue problem on the 
n-dimensional rectangular parallelepiped D: 



(4.12) 



on To, u= §^ = on T'", 




on r- (Fg UP"), 
— X'^Q^u = on Fg, Q — constant > on Tg. 



Similarly, (j4.12p has the special solution u = (^X{xi,--- , .t„_i)) Z(a;„) with Z{xn) 
having form (14. 7p . According to the boundary conditions of (|4.12l) . we get that the 
problem (|4.12p has the solutions 



u{x) = u{xi,- ■■ ,Xn) 



miTT 

C I COS — ; Xl 



COS—; Xn~l Z{Xn), 



ln-1 
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where mi, 
(4.13) 



, m„_i are whole numbers, and Z{xn) is given by 
Z(.„) = cosh/3x„ r(sinh/3/„)(cosh/?;„) + /3/„ 



+ 



sinh^ /3ln - P^ll 
/3(sinh;3/„)(cosh/3Z„) + P^l 



sinh 



sinh^ /3Z„ - 



Xn cosh /SXn 



/? sinh^ 



sinh^ ;3/„ - PHI 



Xn sinh /3a;„ , 



/? = with Y.Vi ^ 0. Since Z"'(0) = 2/3^ ( ^^'""f "J^;°f/;.l^^'" ) 



and^ 



a;„=0 



- •• ,a;„_i))Z"'(0), we get - A^g^^ = q on T^, where 



2pHl ( (sinh /3Z„) (cosh ^/„) + /3/„ 



sinh^ pin - PH 



2/2 



5. Asymptotic distribution of eigenvalues on special domains 



5.1. Counting function A{t). 

In order to obtain our asymptotic formula, it is an effective way to investigate the 
distribution of the eigenvalues of the transformation Gk. (respectively, Gx.d, Gj\f) instead 
of the transformations G^*^ (respectively, G^*] , cj^^). It follows from ([31121) , ^AB, 
and p.30p we obtain 

(5.1) ^ik = {l + Xl)-\ fc = 1,2,3,. ••, 

where /ifc denote the fc-th eigenvalue of Gic or G^^d or G^, and is the fc-th eigenvalue 
of G^^ or G^] or g]^'. Since A{t) = J2\k<T have 

(5.2) A(r)= J2 1- 

5.2. Z) is an n-dimensional rectangular parallelepiped and g'*^ = (J'''. 

Let D be an n-dimensional rectangular parallelepiped, g^^ = J*'^ in the whole of D, 
g — constant > on one face of the rectangular parallelepiped and g = onTg — r+, 
i.e., D ^ {x e R"|0 < x^ < h, i = 1, - ■ ■ ,n}, r+ ^ {x e R"|0 < x^ < h when i < 
n, Xn = 0}) and Fq = F'" = {a; e R"|0 < Xi < k when i < n, Xn — In}- Without loss of 
generality, we assume h < /„ for all i < n. 

For the above domain D, except for the K{D) and K'^{D) in Section 3, we introduce 
the linear space of functions 

— du 

K°(D) = {u\u e Lip(D) n H^(D), — = on F. U Fq, u = on F - F+}, 

Ol' ^ 

Clearly, 

(5.3) K°{D) C K{D) C K'^iD), 
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Closing K'-\ K and K"^ respect to the norm ||m|| = -^Z {u, u), we obtain the Hilbert spaces 
/C", K. and /C^, and 

(5.4) /C° c /C c IC^. 
According to Theorem 3.3, we see that the bihnear functional 

(5.5) M-f ^^u.ds 

•^re- 
defines self-adjoint, completely continuous transformations G", G and G"' on /C", K. and 
/C^, respectively (cf. Section 3). Obviously, 

{G'^Ujv) = {Gu,v) for all u,v in /C°, 

{Gu,v) — {G'^u,v) for all u,v in /C, 
from which and applying Theorem 1.4 of | 26l we immediately get 

(5.6) f^l < t-^k < f^t fc = l,2,3,---, 

where and {/ij?} are the eigenvalues of and C*, respectively. Hence 

(5.7) A°{t) < A{t) < A'^(t) for all r, 
where 

(5.8) A"iT)= 1 
and 

(5.9) A'^ir)^ J2 1- 

Ai^>(l+r3)-i 

We shall estimate the asymptotic behavior of A'^{t) and A'^{t). It is easy to ver- 
ify (cf. Theorems 3.6, 3.7) that the eigenfunctions of the transformations G" and C^, 
respectively, satisfy 

A2u = in D, 
^ = on r+. 

ou 
dv 
u) 



(5.10) 



|2i on r'", and u = Au = on T - (r+ U T'"), 
- A^e^u = on r+, constant > on r+. 



and 



A^u = in D, 



lo.iij <; = and Au = on T - (r+ U T'"), 

^^^"^ — K^g^u — on r+, g = constant > on r+. 

As being verified in Section 4, the functions of form 
(5.12) u{x) = c ( sm — ^ — XI j ■ ■ ■ I sm — ^ Xn-i ] y(Xn) 
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are the solutions of the problem ()5.10p . where Y{xn) is given by ()4.8p . Since the functions 
in (|5.12p have derivatives of any order in D, it follows from Proposition 3.9 that they are 
eigenfunctions of the transformation G° with eigenvalues (1 + A^^)^^, where 

^3 ^ 2r]^ll / (sinh?7/„)(cosh77/„) 

1/2 

and ?7 = 2^ (— ) 

.4=1 * 

Note that the restriction of u on r„ 



1,2,3,. 



(5.14) u\ =c sm— — ••• sm— Xn-i 

when mi,-- - ,rn„_i run through all positive integers (see, Section 4), form a com- 
plete system of orthogonal functions in L^(rg). It follows from Proposition 3.5 that 
if TOi, - - - , m„_i run through all positive integers, then the functions (|5.12p form an or- 
thogonal basis of the subspace of /C", spanned by the eigenfunctions of G*^, corresponding 
to positive eigenvalues. That is, when toi, - - - ,m„_i run through all positive integers, 
then (1 -t- A^)~^, where A'^ is given by (|5.13p . runs through all positive eigenvalues of G°. 

Similarly, for the problem (|5.1ip . the eigenfunctions {uk} of the operator on /C^, 
corresponding to non-zero eigenvalues, form an orthogonal basis of the subspace of IC^. 
The non-zero eigenvalues of are /i^ = (1 + k^)"-'^, where is the k-th Steklov 
eigenvalue of (|5.1ip . 

In order to give the upper bound estimate of A'^ (r) , we further introduce the following 
Steklov eigenvalue problem 



A^u = in D 
fi = on 

(5.15) <^ 



#i = on r+, u = #i = on T'" 

TT - -^Ti— - on r - (r+ u r' 



^^^"•^ — g^u = on r+, g = constant > on F^. 
Let 7^ be the fc-th eigenvalue of (|5.15p . By Theorem 3.8, we have 
(5.16) < kI for all fc > 1. 

We define 

(5-17) d = T^, Afir)= 1- 



1 + 7^ 



M{>(l+r3)- 



It follows from ([ETC]) and ([CTT)) that 

(5.18) A'^{t) < A^ir) for ah r. 

We know (cf. Section 4) that the problem (15.151) has the solutions of form 
rKTn\ ( \ ( "^i'^ ^ ( "In-lTT \ 

(5.19) U{x) = C I cos — ^ Xi j ■ ■ ■ I cos — ^ Xn-l j ^yXn), 

where mi, - - - , m„_i are non-negative integers with X]"=i^ Z{xn) is given 

by (|4.13p . This implies that if mi, - - - , m„_i run through all non- negative integers with 
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(5.20) 



3 _ 2/33/3 /(sinh/3;„)(cosh/3?„) + /3;„ 



pni V sinh2/3Z„-/32Z2 



-, 1/2 



runs throughout all eigenvalues of problem (|5.15l) . 

We first compute the asymptotic behavior of A^ir). By (|5.17p . (I5.20p and the ar- 
gument as in p. 44 of [37 or p. 373 of [5] or p. 51-53 of 26 , Af{T) ==thc number of 
(n — l)-tuples {mi, ■ ■ ■ ,to„_i) satisfying the inequality 



(5.21) 



20X_ / (sinh/?/„)(coshffl„)+/3/„ \ ^ ^3 
Q^'ll V siTxh^ pin- 



where toi, • • • ,m„_i are non-negative integers with 'Y^^^i nii ^ 0. By setting 



(5.22) 

we see that 



t{s) = 2s 



3 / (sinhs)(coshs) -I- s 
\ sinh^ s — 



lim i(s)/s 

s— f-t-oo 



We claim that for all s > 1, 

2s^ [-5-3 + 3s sinh^ s + 3(sinh3 s)(coshs) - 3s2(sinh s)(cosh s) - 2s^ cosh^ s] 



t'{s) 



In fact, let 



(sinh s — 



> 0. 



-s^ + 3s sinh^ s + 3(sinh3 s)(coshs) — 3s^ (sinh s) (cosh s) — 25^ cosh^ s. 



Then 



61(1) > 0, and 
e'{s) = 4(cosh2s 

> 4(cosh^ s 

> 4(cosh^ s 

> 4(cosh^ s 
= 4(cosh^ s 



9 9 ^ sinh s 

3 sinh s - 3s^ - 

cosh s 

[3 sinh^ s-3s^ - s^] 



(s" - s^) > for s > 1. 



This implies that 6{s) > for s > 1. Thus, the function t{s) is increasing in [l,-|-oo) 
Denote by s = h{t) the inverse of function t{s) for s > 1. Then 

lim . 1. 

t 2 



Note that, for s > 1, the inequalities t{s) < t is equivalent to s < h{t). Hence ()5.2ip is 
equivalent to 



pin < KgXr) 
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which can be written as 

n-l 

(5.23) ^(m,/;,)' < 



n-l 



T.J =0,1,2,- •• with ^TO,;^0. 



We consider the (n — l)-diniensional ellipsoid 

n-l 



< 



TTln 



Since {t) + 1 just is the number of those [n — l)-dimensional unit cubes of the z-space 
that have corners whose coordinates are non-negative integers in the ellipsoid (see, VI. 
§4 of [S]). Hence {t) + 1 is the sum of the volumes of these cubes. Let V{t) denote 
the volume and T(r) the area of the part of the ellipsoid situated in the positive octant 
Zi> i — 1, ■ ■ ■ ,n — 1. Then 

V{t) <A^{t) + 1< V{t) + {n- l)5r(r), 

where (n — 1)5 is the diagonal length of the unit cube (see, [5] or [H]). Since 



F(r) =I?„_i2-("-i)?i.../„_i 



(n-l) 



1 /3 

by h{t) ^ (I) as t — ^ +00, we get that 



as T — +00. 



Note that 

It follows that 

i.e., 

(5.24) 



T(t) ~ constant • t''' 



Af{i 



hm = c^„_i(^7r)-("-i);i • • • Z„_ie"-i, 



A^{t) 



+ I ^n — 1 n — 1 



(s/T6 7r)("-i) ' 



as T — > +00, 



where |r+| denotes the area of the face r+. 
Next, we consider A'^{t). Similarly, 



(5.25) 

is equivalent to 
i.e.. 



'^V^lj f (sinh rjln ) (cosh -gin ) + rjln " \ ^ 3 
Q^l V sinh^ rjln ~ V^ll ^ 



Trln 



, mi = 1,2,3,- •• 
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Similar to the argument for A-f{T), we find (see also, §4 of (5j) that 

#{(mi, . . . , m„_i)| E (t^)' ^ {^^^)\ = 1, 2, 3, . . . } 



2—1 



n— l^n— 1 



I.e., 
(5.26) 



lim 



(^^)("-i)''" 



as T 



Not that g = OondD- r+, by dSH]), ((5l8| . ([5:24)) and jSJll), we have 



(5.27) 



A(r) 



/ ! ^n — 1 

g as as r ^ +oo. 



(^7r)("-i) 



5.3. A cylinder D whose base is an n-polyhedron of M" ^ having n 1 orthog- 
onal plane surfaces and g*''' = S^'' . 



Lemma 5.1. Let D^-?) = T^^^ x [0,y, j 



1, 2, where Tg 



(1) 



, a;„-i) e 



"'^ I > for I < i < n ~ 1, and — 1}; '^"■'^ ^g^'' is an (n — 1)- dimensional 



.0) 



1,2. 



cube with side length I — maxi<i<„_i^i. Assume that 
Assume also that g is a positive constant on t'^/^ j = 1,2. If I < In, then 

(5.28) liiD^'^)>li{D^'^) for k = 1,2,3,- ■■ , 

where {'^liD^-'^))^ (similar to 7^ of iS.23\) in Theorem 3.8) is the k-th Steklov eigenvalue 
for the domain D'^^^ . 



Proof. Let v']^^ be the fc-th Neumann eigenfunction corresponding to a^^' for the 
{n — l)-dimensional domain v'^p , (j = 1, 2), i.e., 

(5.29) 

Put 



0) 



dv 



L - 



where Z^^\xn) is as in (|4.13p with j3 being replaced by y cej^ . It is easy to verify that 
u'"j^\x) satisfies 



(5.30) 



A^^^') = 

,,0) 



du): 



du 



on F 

on F, 



du\- 



U) 
Q 1 

U) 




{In}, 



on (dri'^) X [o,g. 



du 



I ^^-(7{(i^«))V-i^'=0 onF«. 
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with 




(sinh y a-^' In) cosh \J a)^' /„ + 



It fohows from p. 437-438 of [5] that the fc-th Neumann eigenvalue a^^"* for the domain v'^p 

(2) (21 

is at least as large as the fc-th Neumann eigenvalue a.\. for the domain Tg . Recalling 
that 2s shih'-^T-t^^^ ) increasing when s > 1, we get 

,3 / , \ 3 



{iUd^'^)) >{li{D^'^)) , fc = l,2,3,... 



if Z < In- Here we have used the fact that y a^, > 1 since any Neumann eigenvalue 

for has the form Y^^Zi ('T^)^- other words, if I < In, then the number (r) 
of the eigenvalues less than or equal to a given bound for the domain D'^^^ is at most 
equal to the corresponding number of the eigenvalues for the domain Z?^^-*. □ 

Similarly, we can easily verify that the number [t) of the eigenvalues less than or 
equal to a given bound for the n-dimensional rectangular parallelepiped D is never 
larger than the corresponding number for an n-dimensional rectangular parallelepiped of 
the same height whose base is an {n — l)-dimensional cube and contains the base of D. 



5.4. D is a cylinder and g''^ = 6^''. 

Let D be an open n-dimensional cylinder in M", whose boundary consists of an (n — 
l)-dimensional cylindrical surface and two parallel plane surfaces perpendicular to the 
cylindrical surface. Assume that g^^ = J*'^ in the whole of D, that Tg includes at least 
one of the plane surfaces, which we call F^, and that g is positive constant on r+ and 
vanishes on Tg — r+. We let the plane surface r+ be situated in the plane Xn — and let 
another parallel surface F'" be situated in the plane {x G R"|a;n = In}- We now divide 
the plane Xn = into a net of {n — l)-dimensional cubes, whose faces are parallel to the 
coordinate-planes in Xn = 0. Let Fi, • • • ,Fp be those open cubes in the net, closure of 
which are entirely contained in F+, and let Qp+i, ■ ■ ■ ,Qq he the remaining open cubes, 
whose closure intersect F+. We may let the subdivision into cubes be so fine that, for 
every piece of the boundary of F+ which is contained in one of the closure cubes, the 
direction of the normal varies by less than a given angle whose size will be determined 
later. (This can be accomplished by repeated halving of the side of cube.) We can make 
the side length I of each cube be less than In- Furthermore, let Dj , {j — 1, ■ ■ ■ ,p), he the 
open n-dimensional rectangular parallelepiped with the cube F^ as a base and otherwise 
bounded by the "upper" plane surface F'"* of the cylinder D and planes parallel to the 
coordinate-planes xi = 0, • • • , a;„_i = (cf. [25]). 

We define the linear spaces of functions 

K = {u\ueLip(D)nH^(D), ^=0 on F„ur'", it = on F - F„}, 
' Ol' 

K° ^ {u.j\uj e Lip{Dj) r\H^{Dj), ^ = on F^ U F^", and u = on dDj - F^}, 

= ,P) 
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with the inner products 

{u,v) — Jjj{Au)(Av)dR + Jp Q^uvds for u,v£K, 
{uj,Vj)j = (uj, + [uj,Vj]j ^ {Auj){Avj)dR + /^^ g^ujVj ds for Uj,Vj e K°, 



respectively. Closing K and with respect to the norms ||u|| — ^ {u, u) and — 
y/ {uj,Uj)j, we obtain the Hilbert spaces K. and /C^ (j = 1, • • • respectively. Clearly, 
the bilinear functional 



u^v] — I Q'^ uv ds 



Jr, 

define self-adjoint, completely continuous transformations G and on /C and /C° by 

(5.31) {Gu,v) = [u,v] for u,v in /C, 

(5.32) {G^jUj,Vj) . — [uj,Vj\j for and Vj in /C°, 

respectively. By defining a space 
p 

/C" ^ ^ ®/C° = Ml + . . . + lij- e /Cj"} 

with its inner product 

(5.33) {u°y)^j2{u„v,),, 

we find that the space KP becomes a Hilbert space. If we define the transformation G° 
on KP by 

(5.34) G%" = G>i + • • • + G^up for m° = ui + • • • + in /C", 

we see that G" is a self-adjoint, completely continuous transformation on KP . If we put 

p 

(5.35) [u°y]^Y>,,v,]„ 

i=i 

we find by ((O^ — that 

(5.36) (G%°,w") = for ah m° and u° in /C°. 

We now define a mapping of /C*' into K. Let u*' = ui + • • • + Up, u £ Hj, be an element 
of JC^ and define 

(5.37) u = n%°, 

where u{x) = Uj{x), when x e Dj, and m(x) = 0, when x £ D — U^^^-Dj. Then u £ K, 
and thus (|5.37p defines a transformation II" of /Cj* ® • • • ® /C^ into K.. It is readily seen 
that 

(5.38) [n%°,n%°] = foraU m° and u° in /C". 
and 

(5.39) (G%°,w°) (Gn%°,n%°) foraU u° and w'' in /C°. 
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From ()5.38p and (|5.39|) . we find by applying Corollary 1.4.1 of ^ that 

m" < Mfe for k = 1,2,3,- ■■ 

Tlierefore 

(5.40) A°{t) < A{t). 
The definition of G° implies that 

(5.41) G%°C/C^", {j = l,- ■■,?), 
and 

(5.42) G%° = G°u°, when u° € K.°. 

From ([OS]) . ([OTt . ([O^ and Theorem 1.6 of PS^, we obtain 

(5.43) ^°(r) = ^AO(r), 

where A'j{T) is the number of eigenvalues of the transformation on /C^ which are 
greater or equal to (1+t'^)~^. Because Dj, {j = 1, ■ ■ ■ ,p), is an n-dimensional rectangular 
parallelepiped we find by (15.26^ that 



(5.44) A^Jt) ^ L^„_i(^7r)-("-i)|rj|e"-V"-i as r ^ +oo, 



where jFj l denotes the area of the face Tj of Dj. By ()5.43p and ()5.44p we infer that 

p 

(5.45) A"{t) - a;„_i(v^7r)^("-i) ^ |Fj|e"-V"-i as r ^ +oo. 

Next, we shall calculate the upper estimate of A{t). Let Pj,{j = p + 1, ■ ■ ■ ,q), be 
the n-dimensional rectangular parallelepiped with the cube Qj as a base and otherwise 
bounded by the "upper" plane surface F'" of the cylinder D and planes parallel to 
the coordinate-planes Xi = 0, • • • ,x„_i = 0. The intersection Pj n Z) is a cylinder 
Dj, (j = P + 1, • • • ,q), with Tj := Qj n F^ as a basis. Then 



(5.46) D = Y,D, 



g 



We first define the linear spaces of functions 

K'^ = {u\u e Lip{D) n H^{D), 1^ on Fg, u = = on F'")}, 
Kf = {u\u e LtpiD) n H^D,), |H = on F„ u = = on F^"}, {j = !,■■■ ,q) 
with the inner products 

(5.47) {u,v} = J^{Au){Av)dx + J^^g^uvds, 

and 



(5.48) {uj,Vj)j=l {Auj){Avj)dx+ I g^UjVjds, 

r,- 
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respectively. Closing K'^ and K'^ with respect to the norms ||u|| ~ ^ {u, u) and \\uj\\j = 
^y{uj,Uj)j, we get Hilbert spaces /C* and — 1, • • • ,q), and then we define the 

Hilbert space 

(5.49) = {u'']u'^ = ui + ■ ■ ■ + uq, uj e Kj] 

with its inner product 

(5.50) 

The bilinear functional 

define a self-adjoint, completely continuous transformation G'j on /C^ given by 

(5.52) {G'^Uj,Vj)j — [uj,Vj]j for all uj and vj in K-j. 

The self-adjoint, completely continuous transformation on /C'' is defined by 

(5.53) = Gfui + ■■■ + G'^gUg for u'^ = m + ■ ■ ■ + Ug in IC^. 
With 

(5.54) [u'',v^]=J2[uj,v,]„ 
it follows from (jCTIl) . ([53^ — that 

(5.55) (G'^w'^,?;'^) = [u^w''] for all w'^ and v"^ in /C''. 
Now we define a mapping 11 of /C into JC^. Let ?i e K{D), and put 

w'* = Hu = Ul + ■ ■ ■ + Uq, 

where Uj(x) — u{x), when a; G Dj. It can be easily verified that 

(5.56) (nu,nw) = for all u and w in /C. 
and 

(5.57) (Gu,w) = (G'^nu,nw) for aU u and v in /C. 
Combining (I5.56p . (|5.57p and using Corollary 1.4.1 of [26j, we obtain 

Mfe < for fc = 1,2,3, •• • , 

and hence 

(5.58) A{t) < A-^ir). 

From G'^/C^ C /C^^, (j = 1, • • • , q), and G'^u'^ = G^w'' when u'^ e /C^^, we get 

(5.59) A'^(r)=^A^^(r), 
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where A'j{T) is the number of eigenvalues of the transformation Gj on K.^ which are 

greater than or equal to (1 + t^)^^. Also, we define (t) similar to (|5.15p and (|5.17p . 
i.e.. 



M^>(l+r3)-i 

where 7^ is the fc-th Steklov eigenvalue of the following problem 



( A'^Uj =0 in Dj, 

^ = Oonr„ ^^. = ^ = on r;.", 
p^^oiAu^^Q on ai?,-(r,ur;"), 

^^^"^^ — ^'^ Q^Uj = on Fj, Q = constant > on r+. 
From Theorem 3.8, it follows that 

ll < Kfe for all fc > 1, 

and hence 

(5.60) ^^(t) < Ajir) for all t and j = 1, • • • ,p, 

where is the k-th eigenvalue of the transformation Cj. Since Dj, (j = 1, • • • is an 
71-dimensional rectangular parallelepiped, we find from (j5.24p that 

(5.61) 4(T)^u;„_i(^^)-("-i)|F,|e"-ir"-i, (j = l,...,p). 

It remains to estimate A^(T),{j > p + 1). According to the argument in p. 438-440 of 
[5], each of the (n — l)-dimensional domains Tj is bounded either by n — 1 orthogonal 
plane surfaces of the partition (the diameter of the intersection of any two plane surfaces 
lies between I and 3/), and an (n — 2)-dimensional surface of the boundary dTg (see, in 
two dimensional case, Figure 5 of p. 439 of 0), or by 2n — 3 orthogonal plane surfaces 
of the partition (the diameter of the intersection of any two plane surfaces lies between 
I and 3/), and a surface of the boundary dTg (see, in two dimensional case, Figure 6 of 
p. 439 of [5]). The number q—p is evidently smaller than a constant C//"~^, where C is 
independent of I and depends essentially on the area of the boundary dTg. Now, we take 
any point on the boundary surface of Tj and take the tangent plane through it. This 
tangent plane together with the plane parts of dTj bounds an n-polyhedron of E"^^ with 
a vertex at which n—1 orthogonal plane surfaces meet (see. Figure 5 of p. 439 of [S] in two 
dimensions), e.g., if 'd is sufficiently small it forms an (n — l)-dimensional n-polyhedron of 
M" with a vertex having n—1 orthogonal plane surfaces (the diameter of the intersection 
of any two plane surfaces is also smaller than Al), or else an (n— l)-dimensional 2{n— 1)- 
polyhedron of K"~^ (see. Figure 6 of p. 439 of [5] in two dimensional case), the diameter 
of the intersection of any two plane surfaces (except for the top inclined plane surface) of 
the 2(n — l)-polyhedron is also smaller than Al; The shape of the result domain depends 
on the type to which f ^ belongs. We shall denote the result domains by S'j. The domain 
Tj can always be deformed into the domain Sj by a transformation of the form (|2.5I) , as 
defined in Section 2. In the case of domains of the first type, let the intersection point 
of n — 1 orthogonal plane surfaces be the pole of a system of pole coordinates r, 6*1, 
^2, • • • , and let r = /(6'i, 6*2, • • • , On-2) be the equation of the boundary surface of 
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Tg, r = h{0i, 02, ■ ■ ■ , 071-2) the equation of the inchned plane surface of the rt-polyhedron 
of M"^^ having a vertex of n — 1 orthogonal plane surfaces. Then the equations 

f)' - a f)' - ft ft' -ft ^' - ■ ■ ■ ,6'n-2) 

^\ - til, ^2-^2, ■■■ , tln-2 - ^n-2, T - r - 

J{ai,02, ■ ■ ■ ,fri-2j 

represents a transformation of the domain Tj into the n-polyhedron S'^ of M"^^. For 
a domain of the second type, let Xn-i = h{xi, • • • , Xn-2) be the equation of top plane 
surface of the 2(n — l)-polyhedron and let Xn-i ~ fi^i, ■ ■ ■ ,a^n-2) be the equation of 
the boundary surface of Tg. We then consider the transformation 

, _ , _ , _ h{xi, ■■■ , Xn-2) 

Xi — Xi, ■ ■ ■ , 2; 2 — Xn — 2, ■^n — 1 — ^n — 1 r/ \ ■ 

j(Xi, ■ ■ ■ ,Xn-2j 

If we assume that the side length I of cube in the partition is sufficiently small, and 
therefore the rotation of the normal on the boundary surface is taken sufficiently small, 
then the transformations considered here evidently have precise the form (j2.5|) . and the 
quantity denoted by e in (j2.5p is arbitrarily small. From Corollary to Theorem 10 of 
p. 423 of "P, we know that there exists a number S > depending on e and approaching 
zero with e, such that 

akiS'j) 
ak{Tj) 

where ak{Tj) and akiS'j) are the fc-th Neumann eigenvalues of Tj and S'j, respectively. 
According to the argument as in the proof of Lemma 5.1, we see that 

UiE,))' = -^i(/„a,(r,)), UiE'^))' = t(/„a,(5j)), 

where t{s) is given by ([O^ . and {-f[{Ejjf and {li{E'j)f (similar to 7^ of 

are the k-i\i Steklov eigenvalue for the rt-dimensional domains Ej — Tj x [0, /„] and 

E'j — S'j X [0,^„], respectively. Recalling that the function t = t{s) is continuous and 

increasing for s > 1, and lims-s.+oo — \t S^t that there exists a constant 6' > Q 
depending on e approaching zero with e, such that 



1 



< 5 uniformly for all k, 



- 1 



<,5'. 



In other words, the corresponding fc-th eigenvalues for the n-dimensional domains Ej = 
Tj X [0, In] and E'j = S'j x [0, /„] differ only by a factor which itself differs by a small 
amount from 1, uniformly for all k. Therefore, the same is true also for the corresponding 
numbers A-^ . (r) and AL, (t) of the eigenvalues less or equal to the bound . 

The domain E'j is either a cylinder whose base is an n-polyhedron of R"~^ having 
(n — 1) orthogonal plane surfaces with its largest side length small than Al or a cylinder 
whose base is a combination of such an n-polyhedron of R"~^ and an (n— l)-dimensional 
cube with side-length smaller than 3^; it follows that if I is taken sufficiently small, the 
number A;^ .(t) from some t on satisfies the inequality 

^£;,(^) < Cir-V"-i + C2r-V"~2 

where Ci, C2 are constants, to be chosen suitably. Thus, (r) can be written as 
A{^,{t) = 6'(C3r-V"-i + C4Z""V"-2), where denotes a number between -1 and +1 
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and C3,C4 are constants independent of l,j and r. It follows that 

^ 4^. (r) = r"-l [eC,{q - p^-' + eC,{q - p^-' 1] . 

As pointed out before, {q — < C; therefore, for sufficiently small I, (q — is 

arbitrarily small and we have the asymptotic relation 

(5.62) lim ^^-^(0. 

where ^(Z) — as / — > 0. For, we may choose the quantity I arbitrarily, and by taking a 
sufficiently small fixed I, make the factor of r"^^ in the above equalities arbitrarily close 
to zero for sufficiently large r. Since 

(5.63) (r) < 4 (r) for j = p + 1, • • • , q, 
we get 

(5.64) hm — < hm — =^{l). 

r— )-+oo T r— !-+oo t 

Combining (lOO)) . (jOS)) . ([STSH)) . (|5J0)) . ((SleTt . ((OS)) . (jOSl and (jOi)) we obtain 

(5.65) c.„_i(^^)-("-i)^,"-if]|r,| < hm 44 <Th^4rT 

< L„_i(^7r)-("-i)e"-if]|r,|) + ,(/). 

Letting Z — ?> 0, we immediately see that 'YTj=i 1^^! tends to the area \Tg\ of Tg and 
lim/^o "^(0 = 0- Therefore, (|5.65l) gives 

(5.66) A{r) ^ 3^"-; |r,|g-V"-^ as r ^ +cx), 

(V iDTr)'."^-^-' 

or 

(5.67) A{t) -^^^i- / q'^-Hs as r ^ +oo. 

In the above argument, we first made the assumption that the boundary dT g of T g 
is smooth. However, the corresponding discussion and result remain essentially valid if 
dV g is composed of a finite number of {n — 2) dimensional smooth surfaces. 



6. Proofs of main results 

Lemma 6.1. Let g*' and y"' be two metric tensors on manifold M. such that 
(6.1) |/-/'|<e, = 
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and 
(6.2) 



1 d 



\9\f' 



1 d 



< e, i,l — 1, - ■ ■ ,n 



for all points in D, where D is a bounded domain in M (see, Section 3). Let 

A*i > > • • • > A*ri > • • • > and /i'l > > • • • > /iJi > • • • > 
he positive eigenvalues of G and G' , respectively, where G and G' are given by 

{Gu,v) — / g'^uvds, for u and v in K., 



{G'u,v)' = / uv ds' , for u and v in JC' . 



Then, for fc = 1, 2, 3, 
(6.3) 



c{(l + eM),(l + eM)i/2}) 
(1 + eMY+ 



< 



(min{(l ~ cM), (1 + eM)-i/2})'"=' 
where M and M are constants depending only on g, g' , dg, dg' and D. 



Proof. It follows from ()6.ip that there exists a positive constant M independent of e 
and depending only on g*-' , g'^^ and D such that 

n n n 

(1 + eM)-i 9%ti < E a'^'Uti < (1 + eM) ^ g'h,ti 

■i,/— 1 i,/— 1 i,/— 1 

for all points in D and all real numbers ^i, • • • , in- Thus we have 
(1 + 6M)-"/2/M < < (1 + eMy/^^l 

which implies (see, p. 64-65 of [26]) that 
(6.4) 
and 



(6.5) 
Thus 
(6.6) 
Putting 



(1 + eM}-"/^dR < dR' < (1 + eM)"/^dR 
(1 + eM)-("+i'/2^s < ds' < (1 + eM)("+^)/2^s. 
(l + eM)~("+^)/^[u,w] < [u,u]' < (1 + eAf)("+i)/2[y^u] 



^ii - 9 - 9 1 



%i = ^l^{V¥\9"') 



we immediately see that 



max|w,;;| < e and max|0i;| < e. 
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; + 



Thus, for any u e /C(-D) or u G JC'^{D), we have 

i,l=l * ' 

so that 



E 



E 



dxi dxi 



dxi 



du 
dxi 



9\9 



du 
dxi' 



It follows that 



\Ag,u - Agu\ < e{Mi\Wlu\ + M2\WgU\), 



where |VgU| is defined in an invariant ways as 



dx^dx^ 



du \ f d'^u 



du 



dx"^ J \dxPdxi P'^dx^'J' 



and Ml and M2 are constants depending only on g, g' , dg, dg' and D. Thus, 



(6.7) 

Set 
(6.8) 

(6.9) 



J \Ag,u - Agu\^dR < 2e^ (^M^ J jV^updi? 

Id |A,zf di? 



\Wgu\^dR 



A?(D) = 



inf 



inf 



In \^9v\'dR 



and the spaces K{D) and K'^{D) are as in Section 3. Furthermore, set 



(6.10) 



(6.11) 



ef{D) = 



inf 



Id Wv\'dR 

vGK{D)^ f2\vlv\^dR=l Ijj \Vlv\^dR' 

Id \^9v?dR 

vf^K-'iD)^ f^\Vlv\^dR=l J^\Vlv\^dR' 



Clearly, A°^{D) > Af{D), and e?(L') > QfiD). Similar to the proofs of Lemmas 2.1, 2.2, 
it is easy to prove that the existence of the minimizcrs to (|6.9I) and (|6.1ip . respectively. 
Therefore, we have that Af{D) > and Qf{D) > (Suppose by contradiction that 
Af{D) — and Qf{D) = 0. Then AgU = in _D for the corresponding minimizer 
u G K'^{D) in two cases. By applying Holmgren's uniqueness theorem for the minimizer 
u E K'^{D) in each case, we immediately see that u = in D. This contradicts the 
assumption J^j jVgUpdi? = 1 or jVgUpdi? — 1 for the minimizer u E K'^{D) in the 
corresponding cases). Combining these inequalities, we obtain 



and 



\Ag,u- Agu\'^dR < 2e^ 



\Ag,u - AgU^^dR < 2e^ 



Ml 



Mi 



Ml 
QfiD) ' A 



D 



\AgU\ 



for ueK{D) 



/ |A„uP for u e K'^{D). 
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Thus we have that, for all u G K{D) or u e K'^{D), 

(1-eM) / \Agu\^dR< [ \Ag,u\^dR< {1 + eM) [ \Agu\^dR, 
Jd Jd Jd 

where M is a constant depending only g, g' , dg, dg' and D. Combining this and (|6.4|) 
we get 

(l-eM)(l + eM)-5 / \Agu\^dR< [ \Ag,u\^dR' < {1 + eM){l + eM)'i [ \Agu\^dR. 

Jd Jd Jd 

That is, 

(6.12) {1 ~ eM){l + eM)-^ {u,uy < (it,u)'^ < (1 + eM)(l + eM)* (u, m)^. 
By (inSl) and (|ST^ we obtain that, for ah u e K{D) or m e K'^{D), 

(l + eM)-("+i)/2[„^u] ^ [u,uy 



( max{(l + eA/)(l + eM)t , (1 + eM)("+i)/2}) ((u, u)* + [u, u]) {u, u)'* + [u, u]' 

^ (1 + £M)("+i)/^[u,m] 

~ (niin{(l - eM)(l + eAf)-t , (1 + eM)-("+i)/2}) ((u, m)* + [u, u]) ' 

which implies (|6.3p . □ 

Remark 6.2. Let F and F be two bounded domains in R"^^, and let F is similar to F 
(in the elementary sense of the term; the length of any line in F is to the corresponding 
length in F as ft- to 1), and let Fo = F x {a} and Fq = F x {ha}. It is easy to verify that 

Af{D) = h-^Ai{D), ef{D) = ei{D), 

where D = T x [Q,a], D = T x [0, ha], and Af (D) and &i{D) are defined as in (HH) and 
(|6.1ip . respectively. 

Lemma 6.3. Let G and G' be the continuous linear transformations defined by 
{Gu,v) = / Q^uvds for u and v in K{D) or K'^{D) 



{G'u,v)' = [ g'^uvds for u and v m K{D) or K'^iD), 
Jr, 

respectively. Let 

Ml ^ M2 > • ■ ■ > MA; > ■ ■ • > and Mi 5^ M2 ^ ' ' ' ^ Mfe ^ ' ' ' > 

be the positive eigenvalues of G and G' , respectively. If g < g' , then 
(6.13) Mfc<Mfc /or fc = l,2,3,--- . 



Proof. Since g < g', we see that for any u e K{D) or K'^{D), 

{Gu, u) _ /r, g^"^ ^ /r, ^'^^^ _ {G'u.u)' 

{u,u) {u,u)* + Jp g^u^ ds ~ {u,u)* + Jp g'^u'^ ds {u,u)' ' 

which implies (|6.13p . □ 
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Proof of Theorem 1.1. a) First, let {A4,g) be a real analytic Riemannian manifold, 
and let the boundary dD of D be C^'^-smooth. As in [26], we divide the domain D 
into subdomains in the following manner. It is clear that dD is the union of the por- 
tions Fi, • • • , Fp (without common inner point on the surface). Let U he a coordinate 
neighborhood which contains Tj, let Xi = Xi{Q) and ai = ai{vQ) be the coordinates of a 
point Q in T j and the interior Riemannian normal vq at Q, respectively. We define the 
subdomain Dj and surface FJ by 

Dj = {P\x{P) = x{Q)+S,na{uQ), Q e F,, < e„ < a} 

and 

FJ = {P\x{P) = x{Q) + a a[vQ), Q e F,}, 
where cr is a positive constant. The closure of Dj is 

(6.14) Dj = {P\x{P) = x{Q) + ^naiiyq), Q € f,, < < a}. 

By the assumption, each Tj, which is contained in a coordinate neighborhood, can be 
represented by equations 

(6.15) Xi = V'i(6, • • • ,Cn-l) 

with C^''^-smooth functions ipi^ i-C-i it is the imagine of the closure Tj of an open domain 
Tj of R"~^. Hence, if a is sufficiently small, the definitions have a sense and the formula 

(6.16) x{P)^x{Q)+^na{i^Q), QeTj, 0<^n<<J 

defines a C^'^-smooth homeomorphism of a neighborhood of the image of Dj in M" given 
by the coordinates x and a neighborhood Uj of the closed cylinder Fj in E" defined by 
Fj — ■ ■ ■ ,^n-i) G Tj, < ^„ < a}. (Therefore all Lemmas of Section 2 are true 

for every Dj). Moreover, the domains Z)i, • • • , Dp have no common inner points and the 
remainder Dq — D — U^^^^Dj of D has a finite number of connected parts. Note that the 
boundary of Dq contains no part of Tg with measure > 0. 

Define the space N — N{D), Af and the transformation G on A/" as in Section 3. We 
shall investigate the asymptotic behavior of A{t) with regard to transformation G on 
space Af. Moreover, we define the function spaces 

= {uj\uj e Lip{Dj) n H'^{Dj), ^ = on F^- U FJ, and u = on dDj - Tj}, 

HS = {uo\Lip{Do) n H^Do), u = ^=0 on dDo}, 
Kf = {uj\uj e Ltp{Dj) n H^Dj), ^ = on F„ u, ^ ^ = on FJ}, 

(i = 0,1,--- ,p), 

and the bilinear functionals 

(6.17) {u,,u,))= [ |AgM,fdi?, (j = 0,l,--- ,p), 

J Dj 

(6.18) [uj,Vj]j^ g^UjVjds, (j = l,---,p), [uo,vo]=0, 
and 



(6.19) 



{uj^Vj)j = {uj,Vj)*j + [uj,Vj]j, {j = 0, 1, • • • ,p). 
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where Uj,Vj £ or K!j. Closing and Kj with respect to the norm \uj\j = 
■y/ {uj, Uj)j, we get the Hilbert spaces K° and /C^^, {j 0, 1, • • • Then, in the same 
manner as in Section 5 we can define the Hilbert KP and JC^, and define the positive, 
completely continuous transformations G"^, C^, and Gj on /C°, /C^, JC'j and /C^, re- 
spectively. Furthermore, we can prove 

(6.20) A°(t) < A{t) < A''{t) for all r, 
and 

(6.21) A°(r) = A°(r), A'^(r) = A^^(r), 

where ^"(t), A'^(t), ^^(t) and Aj{t) are the numbers of eigenvalues of the transfor- 
mations G^, C^, G° and Gj on /C'', /C'', /C^ and ICj which are greater than or equal to 
{1 + T^)~^, respectively. 

Since [uo,uo]o = for all UQ,va e or K^j and (GS]uo,uo)o = (GoWo,mo)o = [wo,-uo]o, 
we find fairly easily that Gg = GJ^j = 0, so that A'^ir) = A[5(t) = 0, (t > 0). Thus, we 
have to estimate A^{t) and ^^(t) for those domains Dj, where g^ds > 0. 

We can choose a finer subdivision of dD by subdividing the domains Tj into smaller 
ones, e.g. by means of a cubical net in the coordinates ^. According to p. 71 of [55], by a 
linear transformation of the coordinates we can choose a new coordinate system r] such 
that 

/(^)=5'', (z,Z = l,--- ,n), 
for one point in the mapping Tj of Tj , and 

(6.22) a;fc(P) = ipkivir ■ ■ + Vna{v{r]i, ■ ■ ■ ,?7„_i)), 

for (771, •• • , Vn-i) £Tj, <r]n <cr 

defines a G^'^-smooth homeomorphism from Ej to Dj, where Ej = {rj = {rji,- ■ ■ , ?7„)| (t?i, 
• • • , rjn-i) € Tj, < 77„ < cr} is a cylinder in M". (This can also be realized by choosing 
a (Riemannian) normal coordinates system at the point fj e Tj for the manifold {Ai,g) 
(see, for example, p. 77 of [H]) such that a{v{ri)) = (0, • • • , 0, 1) and by using the mapping 
(16.221) ). If we denote the new subdomains of dD by Tj as before, it is clear that we can 
always choose them and a (see p. 71 of [H]), so that, 

(6.23) \9"\^') ~ 9'\n)\ < ^, z,; = i, 



3.24) 



1 d , „, 1 d 



i,l — 1, - ■ ■ ,n, 

for any given e > 0, and all points rf e Ej. The inequalities (|6.23p imply that 



n 

2 



3.25) (1 + M,e)-i Y.*^<J2 9"iv%ti < (1 + M,e) ^ t 



i = l i.l = l 1=1 



for all points 77' £ Ej and all real numbers ti, - ■ ■ ,tn, where Mj is a positive constant 
depending only on 5'' and Ej (cf. Lemma 6.1). This and formula (128) of [2^ say that 

(6.26) (1 + 6A-f,)-"/2|T,| < |r,| < (1 + eM,r/'\T,\, 
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where 

\^j\= Vaiv) dr]i---dr]n-i, \Tj\ = dr]---dr]„_i 

are the Riemannian and Euclidean areas of Tj and \Tj\, respectively, 

Now, wc consider the Hilbert spaces /C° and JCj. When transported to Ej, the under- 
lying incomplete function spaces Kj and are 

K° = {u\u e Lip{Ej) n H^{Ej), ^ = on TjU T/, and u = on dEj - Tj} 
and 

Kf = {uj\uj e Lip{Ej) n H'^iEj), = on Tj, uj = = on T/}, 
respectively. The inner product, which is similar to Section 5, is defined by 



{u,v)j^ / {Agu){Agv)y^g{T])drji ■ ■ ■ drj^ + / q^uv^J g{ri) drji- ■ ■ rjn-i 
and the transformations G° and Gj are defined by 

{G%,v)j= / g^uv^/g(r])dr]i---rin-i, for u,v in K.j, 



and 



{GjU,v)j = / B^uv\/g{ri) drji-- ■ r?„-i, for -«,i> in IC 



respectively. 
Let 

(6.27) Q . = inf g and gj — sup g, 

and let us introduce the inner products 

{u,v) . = / (Au)(Ai;)rfr/i • • • d?7„ + / ^.uv drji ■ ■ ■ drjn-i 

and 



(m, w)^- = / {Au){Av)dT]i ■ ■ ■ drjn + / eifui; drji--- drjn-i 
in the spaces K^^ and K"^, respectively. By closing these spaces in the corresponding 



norms, wc get Hilbert spaces ]C_j and JC^ . Furthermore, wc obtain the positive, completely 
continuous transformations and on K^^ and /C^, which are given by 

(6.28) {G!jU,v)j = / g^.uv drji ■ ■ ■ drjn-i, for u and v in £^ 
and 

(6.29) {G^u,v)j = / gjuvdrji- ■ ■ dr]n-i, for u and i> in 
respectively. 
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Let /ifc(G°) be the k-th positive eigenvalue of G° and so on. According to Lemma 6.1 
and Remark 6.2, Af (Dj) and Qf{Dj) have uniformly positive lower bound when repeated 
taking finer division of D (In fact, by repeated halving the side length of every rectangular 
parallelepiped in the partition net of the coordinates rj for each cylinder Ej , we see that 
Af (Dj) will tend to +oo, and that 'df{Dj) will have a positive lower bound). This implies 
that the corresponding positive constants Mj and Mj have uniformly upper bound when 
we further divide the domain D into finer a division, where Mj is defined as before, and 
Mj is a constant independent of e and depending only on g, dg and Ej as in Lemma 6.1. 
Denote by Cj{e) the maximum value of (1 + eiVf,)"'^' ( min{(l — eM), (1 + eM)^^}) ^ 
and (1 + eMj)"+^(max{(l + eM), (1 + eM)i}). Obviously, Cj{e) 1 as e ^ 0. By 
virtue of (|6.23|) and ()6.2|) . it follows from Lemmas 6.1 and 6.3 that 

(6.30) MG'^) < cj{e) ^ik{G'^) 
and 

(6.31) Mfc(G°)>c,(e)-Vfc(G°), 
so that 

A'^iT)<A'^{c,{e)T + c,{e)-l) 

and 

(6.32) A,"(r) > 4(c,(6)-V + c,(e)-i - l) 

where (t) and iij'(T) are the numbers of eigenvalues of the transformation Gj and G° 
which are greater than or equal to (1 + t^)^'^, respectively. By ()6.20|) and ()6.21|) . we 
obtain 

(6.33) ^4(c^.(e)-V + c,(6)-i-l) < A(r) 

j 

j 

Finally, we shall apply the results of Section 5 to estimate i4°(r) and A'^{t). Note 
that 

(6.34) 3J(t) < 1J(t) for all r > 0, 

where is defined similar to ((??T5]) — ([CTT)) . Formula dOI]) . !i^M\ and imply 
that 

(6.35) hm M^>c.„_i(^^)-("-i)|r,|e"-i 
and 

(6.36) li^ <^„_,(^^)-("-i)|r,.|^-i, 

where \Tj\ is the Euclidean area of T,-. Combining (H^Ol), 
(jOTI) . ((04)) and ((5:65)) we find 

Ih^^(T)T-("-i) <c.„_i(^7r)-("-i)5,(e) Ve"-i|r,|, 
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and 

lim A(T)r-("-i) >c^„_i(^7r)-("-i)£,(e)-i^e';-i|r,|, 

r — >oo . 

3 

where Cj(e) = (1 + eA/j)"/^Cj(e)"^^. Note that g is Riemannian integrable since it is 
non-negative bounded measurable function on dD. Therefore, letting e — > 0, we obtain 
the desired result that 

(6.37) ^(-)-(^ff^/ ^""^^ as.^+oo. 



dD 

b) Since a metric g*' in can be approximated by a metric g'^' which is C'^-smooth 
on M and piecewise C^^'^-smooth in any compact submanifold of {Ai,g) such that 

|gf -/!<£, i,l = l,---,n, 




< e, 



i,Z = 1, 



,n, 



for all points in D, with any given e > 0. In addition, any bounded domain D with C^- 
smooth boundary can also be approximated (see, the definition in Section 2) by domain 
D'^ with C^-smooth and piecewise C^-'^-smooth boundary. Thus, by Lemma 6.1 and a), 
we can estimate the eigenvalues for g'J-^ in D'^. But for these eigenvalues (|6.37l) is true, 
so that letting e — and noticing that ds'^ — > ds, we get that (|6.37p also holds for such 
metric g'' and D. □ 



Proof of Corollary 1.2. By (jl.Sp . we have 

Wri-lAj, 



(6.38) 



vol(9I?)), as fc — >■ +00. 



(^7r)(«-i) 

Since A{Xk) — k + 2, we obtain (|1.6p . which completes the proof. □ 

Proof of Theorem 1.3. Let T, : H^/^{dD) -> H-'^/^{dD) be defined as follows: For 



any cj) e H^/\dD), we put T,(^ ^ [j^ 



£)^ diy 



\dD 



where u satisfies 



A^u = in D, 
§H = on dD, 
u - 



on dD, 



and e > is a sufficiently small constant. Clearly, is a self-adjoint, elliptic, pseudo- 
differential operator of order 1. Also, its principal symbol has the form c{q + e)~'^|C|^, 
where ^ € M"~^ and c is a unknown constant (we will determine it later). It is easily seen 
that the operator has the same eigenvalues (Afe(e))'^ (respectively, the corresponding 
eigenfunctions uu) as the following Steklov eigenvalue problem: 



Aluk = in D, 
^ = on dD, 

%^-(A.(6))'(^.- 



e)^Mfc = on dD. 
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Applying the standard method (see, for example, p. 163 of [27], [H], [32] or [2]), we have 



(27r) 



n-l 



an Jo 



(27r)"-i 7a^V ci/3 

where ^'^^(r) = #{(Afe(e))'^ < t'^}, o-„_i and — are respectively the area and 
volume of the unit sphere and the unit ball of R"~^. Letting e — 0, we obtain 



It follows from the leading asymptotic formula (|1.5p of Theorem 1.1 that the constant c 
must be 2. Therefore, 

A{t) = 2-(t^) • (27r)-("-i'a;„_ir"-i / p''-\s)ds + 0(t("-2)) as t +00. 

□ 

Remark 6.4. It is worth noting that for the harmonic Steklov eigenvalue problem 
(|1.2p with smooth boundary dD, we can also obtain a better asymptotic formula with 
remainder estimate than that of [26] (i.e., (|1.3p ). In fact, we first define the pseudodif- 
ferential operator Se : H^/'^{dD) H-^/'^{dD) as follows: For any e H^/'^{dD), we 
put St<j) '■— ((f? + ^)~^^) where v satisfies 

AgV = in £), 



V ~ (j) on dD. 

This is just the well-known Dirichlet-to-Neumann map, and it has the same eigenvalues 
as the Steklov eigenvalue problem (|1.2p with g being replaced by p + e. One knows (see, 
for example, p. 103 of [29]) that the principal symbol of the pseudodifferential operator 
T is (p + e)~^|CI- So, we have the following asymptotic formula (see, for example, p 163 

of [13, m) 

B('\t) = (27r)-("-i) ( / ds I an +0(r("-2)) 

\JdD J (Q+'^)-^\i\<r J 



(27r)"-i 
By letting e 0, we get 



(q{s) + ey-^ds + 0(t("-2)) as t ^ +00. 

dD 



(6.39) B{t) = ^ / e"-i(s)ds + 0(r("-2)) as r ^ +00. 
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